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\^ Abstract 

a^ ■ 

I We establish and discuss a priori estimates on subgrid stress and subgrid flux for filtering 

■ schemes used in the turbulence modelling method of Large- Eddy Simulation (LES). Our 

O \ estimates are derived as rigorous consequences of the exact subgrid stress formulae from 

' Navier-Stokes equations under realistic conditions for inertial-range velocity fields, those 

conjectured in the Parisi-Frisch "multifractal model." The estimates are shown to be an 

' expression of "local energy cascade," i.e. the dominance of local wavevector triads in the 

H \ 

energy transfer. We prove that for nearly any reasonable filter function the LES method 
defines an energy flux in which local triads dominate in individual realizations, due to 
cancellation of distant triadic contributions by detailed conservation. A somewhat similar 
observation of Leslie and Quarini on graded filters in the EDQNM closure is shown to be 
unrelated to the cancellation we establish in Navier-Stokes solutions. The sharp Fourier 
cutoff filter is one example which does not satisfy the modest conditions of our proof and, in 
fact, we show that with that filter the energy transfer in individual realizations at arbitrarily 
high Reynolds number will be dominated by nonlocal, convective sweeping. 

Key words: Navier-Stokes, turbulence, large-eddy simulation, multifractal model 
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1 Introduction 



A fundamental hypothesis in the Kolmogorov phenomenology of developed turbulence is the 
"local energy cascade" 0, ^. In the first of his 1941 papers, Kolmogorov expressed this 
in spatial terms as a process of energy transport to small scales by successive transfer from 
"pulsations" of one size to daughter "pulsations" of size about half as large. The parallel 
works of Heisenberg Q, von Weizsacker and Onsager |^] stated this in spectral terms 
as a process in which energy is transported through Fourier space along chains of wavevector 
triads in which typically two wavevectors of the same size transfer energy to a third wavevector 
about twice as big. Recently, the dominance of these local wavevector triads to turbulent 
energy transfer has been the subject of much debate. The papers P|-|2C] in chronological 
order are a representative sample containing many further references. The cited studies are 
based mostly upon direct numerical simulation of the Navier-Stokes equations or upon the 
solution of analytical closures, such as the EDQNM transfer equation. The first method is 
severely constrained by the low Reynolds numbers attainable with current computer resources. 
The numerical solution of analytical closures, on the other hand, is feasible at high Reynolds 
number — at least in the case of homogeneous and isotropic turbulence — but is subject to doubt 
concerning the uncontrolled closure hypotheses invoked. 

However, we have in fact demonstrated in our past work |^l| locality of instantaneous energy 
transfer for individual turbulent solutions of the Navier-Stokes equations based upon a direct, 
rigorous analysis of the inertial terms and thus valid for Reynolds numbers arbitrarily large 
without any assumption of statistical homogeneity or isotropy. Our original study involved 
mathematical considerations on so-called "weak solutions" which might make it inaccessible to 
fluid-dynamicists without a good grounding in modern analysis. Therefore, we shall present 
here the basic arguments of our earlier work with an emphasis on their physical interpretation 
in order to place them in the context of the existing turbulence literature. In addition, we shall 
discuss stronger results on simultaneous space-scale locality of transfer which we subsequently 
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obtained and new material exposing the various relations between both of our two works 
and the previous literature. 

We have found that the filtering technique commonly used in the Large-Eddy Simulation 
or LES modelling method (see 1 23 , 24 1 ) provides the most theoretically well-founded definition 



of the intuitive idea of "scales of motion" in turbulent flow. In this method a filtering function 
G£{r) = £^'^G£{r/i) is convoluted with the velocity field, as 

v,(r)^ J (iVG,(r-r')v(r'), (1) 

to define the "large-scale velocity" or "resolved field" v^. The "small-scale velocity" or "subgrid 
field" is then simply defined as the complementary component: 

vKr) = v(r)-v,(r). (2) 

When the filtering operation is applied to the Navier-Stokes system an equation is obtained for 

dtVe + V • {vfVi + Te) = -Vpi + 1^0 A vg, (3) 
in which p£ is the filtered pressure field (required to maintain V ■ = 0) and 



Te = (vv)^ - vgvi, (4) 

is a tensor representing the "subscale stress" of the eliminated turbulent eddies. 

We shall use the same filtering technique as in conventional LES but in order to derive 
rigorous mathematical estimates for the exact subgrid stress which follows from Navier- 
Stokes. We show under very mild assumptions on the filter function G that the energy fiux 
which appears in the LES method is dominated by local wavevector triads, in agreement with 
the conventional wisdom. However, a few commonly used filters — including the sharp cutoff 
filter in Fourier space — will be shown to violate these assumptions and to lead to formulae 
for subscale stress which are not appropriate to model as an effect of the small-scales alone. 
The LES technique with a "good"' filter will in fact be shown to simply and naturally resolve 



certain subtleties of the locality issue, discussed in which arise purely as an artefact of 

the conventional use of Fourier series. 

Our results have also some further implications. First, the derived estimates place an a priori 
constraint on subgrid models and filters and may be used to evaluate the various choices. This 
is the subject of a companion paper |25], hereafter referred to as II. Second, the estimates imply 



restrictions on the space-regularity of turbulent velocity fields in the inertial range of scales in 
order to be consistent with the basic phenomenology of constant mean energy flux. Roughly 
speaking, the velocity must appear on those scales as a continuous, nowhere-differentiable field 
of the type represented by the famous Weierstrass function or a typical path of Brownian 
motion. This was noted already by Onsager in 1949 [^. 

In fact, all of the estimates we derive are based upon the regularity of inertial-range velocity 



fields postulated in the Parisi-Frisch "multifractal model" [26|, who generalized Onsager's ob- 
servation. See 1^,1^ for recent accounts. Let us just recall here that those authors proposed 
that in the limit of zero viscosity the velocity field in turbulence should be Holder continuous 
at each point, that is, for every point r in the flow domain an inequality 

|v(r + l)-v(r)| <C|1|'' (5) 

should be satisied for some real h and all |1| < 1. (This definition must be modified for negative 

h.) If so, then to each point may be assigned the largest value of h for which Eq.(|5|) holds and it 

was furthermore postulated that the points where the value < h occurs should be a fractal set 

S{h) with dimension D{h). This is the origin of the word "multifractal." Notice that the original 

1941 Kolmogorov theory corresponds in this language to D(l/3) = 3 and D(h) = —oo otherwise, 

that is, /i = 1/3 everywhere. Q Although the model must be considered to be conjectural, we 
^The original motivation for this "multifractal hypothesis" was the experimental observation that longi- 
tudinal velocity-difference moments (structure functions) scale as powers of the separation length £ — 



SAe) = ([i-(v(r + l)-v(r))]") 

X 



<(tV'. (6) 
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shall see that it has some mathematical support from the Navier-Stokes equations. In fact, the 
"multifractal" velocities are the most regular velocity-fields possible which are still compatible 
with the basic physics of constant mean energy flux through inertial-range scales. In particular, 
it is not possible for the velocity to be everywhere space-differentiable on inertial-range length- 



h > 1/3 at all r and inertial-range £. The argument for this is simple, but somewhat off the 
main line of the paper, so that we include it as an Appendix. 

The contents of our paper are as follows: In the next Section II we shall derive a priori 
estimates on subgrid stress and flux as defined in the LES method, which are exact consequences 
of Navier-Stokes assuming the "multifractal" regularity of velocities. In Section III we explain 
how the derived estimates express the "locality" of the instantaneous energy-transfer. It is 
shown by explicit analysis that for a "good" choice of filter the LES flux gets a negligible 
contribution from distant wavenumber triads and we thereby resolve some of the controversies 
regarding "locality" of transfer. We contrast our result with an earlier observation of Leslie and 
Quarini . We also show that there are "bad" choices of fllter for which the flux is dominated 
by distant triadic interactions of convective type. Finally, in Section IV we make some closing 
remarks and conclusions on the "multifractal model," the locality of transfer, and the filtering 
approach. 

over inertial-range scales <^ £ <^ L and that C^p 7^ p/3, the prediction of K41 theory. Parisi and Frisch explained 
this latter fact heuristically by noting that, with their hypotheses, 



where the second factor measures the fraction of boxes of size £ in the domain of size L in which the exponent h 
occurs. Applying steepest descent to evaluate this integral for ^ <C L they obtained Eq.(@) with 



scales (which would imply that Eq.(|5|) hold at all points r with h > 1.) Just to be clear, 
we emphasize that the conclusion is that Eq.(||) cannot hold for Navier-Stokes solutions with 




(7) 



Cp ^■mf[ph+ {3 - D{h))]. 



(8) 
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2 A Priori Estimates and the "Multifractal Model" 

We have already introduced the concept of "subscale stress" T£ in the LES filtering technique. 
It is also natural to write down a local energy balance for the resolved kinetic energy 

e,{r,t) = ^vj{r,t), (9) 



as 



dteeir, t) + V • Uir, t) = -n^(r, t) - e^(r, t). (10) 



For example, see Section 6.1 of |24|. Each of the terms in the microscopic balance relation 
Eq.([To|) has a precise physical interpretation. The term e^(r, t) = z^o(Vv£)^ represents the local 
dissipation of energy in the length-scales > £ through the action of molecular viscosity. The 
term j^(r, i) = (p^ + etyv^ + ti ■ V£ — i^oVei is a spatial current of energy in the scales > £. 
Finally, the term of main interest for us 

n,(r) = -(Vv,) : re (11) 

represents a local energy flux from the large-scales to the small-scales. It gives the "effective 
dissipation" of energy in the large-scales > £ due to the action of the eddy-stress of the small- 
scales < i on the gradients of the large-scale motion. 

In this section we shall derive some simple a priori estimates for the quantities T^(r) and 
n£(r) assuming the Holder regularity Eq.(|5|) at point r, a condition usually denoted v S C^(r). 
Only afterward will we consider the consequences of the estimates and the justification for the 
assumed regularity. A key identity which was noted, in a different language, in a paper of 
Constantin et al. Il3l|l, is 



T,(r) = (Av(r)Av(r)), - v;(r)v;(r). (12) 

where 

Asv(r) = v(r) - v(r - s) (13) 
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is the "backward difference" and 

(/), ^Jds G,(s)/(s), (14) 

denotes an "average" over the separation distance with respect to the filter function. The proof 
is elementary. Writing the righthand side of Eq.(^) as 

I ds GKs)[v(r) - v(r-s)][v(r) - v(r-s)] - [v(r) - v,(r)][v(r) - v^r)], (15) 

using the normalization condition 

J d'^v G,(r) = 1 (16) 



and some shifting of integration variables and cancelling terms gives (vv)^ — v^v^, which is the 
definition of the subscale stress. 

A more symmetric and suggestive form of the same identity is 

T,(r) = (Av(r)Av(r)), - (Av(r)),(Av(r)),. (17) 

The proof is instructive. By definition 

v,(r)= j dVG,(r-r')v(r'). (18) 

Therefore, using = v — and the normalization condition Eq.(16) for G^, 

v;(r) = j d^r' Ge{r - r')[v(r) - v(r')] 

= J d'^s G,(s)Asv(r) = (Av(r)),. (19) 

Therefore, the subgrid velocity is represented by an "average" over velocity-differences with 
separation < £. Substituting into the previous formula Eq.(|l2|) for ti gives Eq. (p^ . In this 
expression for the subscale stress the velocity again appears only through its difference Av over 
length-scales < i. If the filter function G is positive as well as normalized, then Eq.(jl^) defines 
a true average and, from Eq. ([T7|), ti is a positive, symmetric tensor. An equivalent identity 
was discovered and used to establish positivity of the stress in a recent work of Vreman et al. 
which appeared after the first submission of this paper. 



We can now straightforwardly derive our a priori estimates using the condition v G C^{r) 
i.e. |Asv(r)| < (const. )|s|''. For example, it follows directly from Eq.(19) that 



|v^(r)| < (const.) J d'^s Gi{s)\s\^ 

= O . (20) 

if the filter function satisfies the modest condition / dx |G(x)| • |x|'^ < oo. Likewise, 

T,(r) = O [f'') (21) 



by using Eq.([T7|) if / dx |G(x)| • jxp'^ < oo. This is our main estimate on the subscale stress. 
To obtain an estimate on the subgrid flux II^ we also require a formula for Vv^ expressing it 
in terms of velocity-differences over length-scales < £. Using / VG = we find 

Vv,(r) = y'dV(VG,)(r-r')[v(r')-v(r)] 

= -l^^s (VG,)(s)Asv(r). (22) 

Making the same estimations as before, we then find that 

Vv,(r) = O (i^-') , (23) 
if / (ix |VG(x)| • Ixl'* < oo. From this result and the formula Hi = — Vv£ : T£ it follows that 

n,(r) = O (f'^-') . (24) 



This is our main estimate for the flux. 

If we gather together all of the conditions on the filter function G, we see that it is sufficient 
for the previous estimates to apply that it obey the moment condition 

y dx [|G(x)| + |VG(x)|] • |x|2 < oo, (25) 

when the Holder indices lie in the range < /i < 1. The condition is quite modest. Of filters 
commonly used in practical LES, the Gaussian filter is a "good" example which easily satisfies 
these constraints. On the other hand, another common filter, the sharp cutoff filter, violates 



condition Eq.(p^), because of slow spatial decay in space ~ |x| for large x. In fact, we will see 



by an example from [21| that the previous bounds for and 11^ may actually be violated for 
this filter. Smoothly cutting off the tails of the filter in physical space would cure this problem 
in principle, but there are definite signs of the bad behavior in practical use of the sharp-cutoff 
filter in current LES simulations. This is discussed at length in II. 

So far we have employed the hypothesis of local regularity made in the Parisi-Frisch "multi- 
fractal model," without considering its validity. Since the assumption of local Holder regularity 
cannot presently be derived a priori as a theorem for solutions of the 3D Euler equations, it does 
require some justification. Mathematical results — presented in — in conjunction with basic 
facts of turbulence phenomenology strongly suggest the "multifractal model" as the correct 
candidate to describe the space variation of inertial-range velocity fields. We must emphasize, 
however, that the key elements of the "multifractal model" used in our previous analysis do 
not depend upon the questionable experimental results on anomalous scaling |^. On the con- 
trary, much weaker conditions guarantee the validity of the local Holder regularity we employed. 
A particular condition which suffices is that the inertial-range energy spectrum be power-law 



E{k) ~ with any spectral exponent whatsoever. In fact, we have proved elsewhere [331 
that any homogeneous random field with such a power-law spectrum has realizations with local 
Holder regularity with probability one. We just emphasize here that "intermittency" is irrele- 
vant for this issue and that the spectral exponent could be the Kolmogorov one, n = 5/3, or 
anything else. Furthermore, local Holder regularity is supported by the estimate Eq . (p4|) , 
which implies H^ — > as £ — > if /i > ^. It is impossible to account for mean energy flux 
through the inertial range if that condition holds (e.g. if the velocities are spatially smooth) at 
all points of the domain on inertial-range scales. Therefore, the "multifractal model" hypothesis 

of local Holder regularity is vital to account for turbulent energy balance in the inertial-range. 

^If the energy spectral exponent is exactly n = | , then "negative Holder singularities" may occur but can 
occupy a space set of Hausdorff dimension equal to at most 2| |^^. If n > |, then the Hausdorff dimension of 
the "negative exponent" set is even smaller. 
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3 Locality of Energy Transfer 

3.1 The Estimates and Local Transfer in Wavenumber 

The a priori estimates on subscale stress and energy flux in LES were derived in the previous 
section just as mathematical upper bounds. However, they have a very intuitive physical 
interpretation in terms of locality of energy transfer, which we discuss in this section. An 
essential point in the LES method is that it defines, via Eq.(|ll|), an energy flux which is 
proportional only to the gradient of the large-scale velocity field. The physics of this expression 
seems transparent: it describes the energy transfer due to the stretching action of the large- 
scale strain = ^ (^Vv^ + Vv^^ upon the stress field T^(r) due to the random distribution 
of vorticity elements in the small-scales < i. In fact, it is often convenient to regard the fine- 
grained turbulence as a "tangle" of vortex filaments with a viscoelastic response to imposed 
large-scale strain |3^. All of the explicit contributions of convective processes to the rate of 
change of the resolved energy ee{r,t) at the point r in the LES method are incorporated into 
the term (e£ + P£)v£ of the current and can be accounted as contributions to transport of 
energy in space rather than "downward" to smaller scales. Since it is the large-scale convective 
processes which can produce apparent violations of locality in energy transfer in the wavenumber 
description, it is rather natural to find that ^e{Tc) obeys a "local estimate" ~ £3h-i ^j^j^j^ arises 
from Vv£ ~ i^-^ and n ~ v^v^ ~ f^. 

However, we shall show that this simple picture actually depends upon the technical-looking 
condition Eq.(|25|) we imposed upon the filter function G in the previous section. Contrary 
to naive expectations, the LES flux can retain implicit contributions from purely convective 
processes if that (mild) condition is violated. We shall show this by a concrete example. In 
such cases it is inappropriate to estimate ri ~ v^v^ because the "subscale stress" is actually 
still dependent upon the large-scale motion! In the course of our analysis we will make contact 
with the traditional discussion of the locality issue in Fourier space and resolve the subtle points 
which arise in that description. 
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In the Fourier representation, the conventional definition of "(subscale) energy flux" n(A;, t) 
goes back at least to the work of Heisenberg [Q] and von Weizsacker who defined it as the 
total output power from the Fourier modes in the sphere of radius k in wavenumber space: 

dE<''{t) 



U{k,t) 



(26) 

Euler 



dt 

Here E^^{t) is the instantaneous energy in modes of wavenumber < k. If the sharp cutoff filter 
is used in LES in a rotationally symmetric form, as 

f 1 if Ifcl < 27r 
G(k) = , (27) 

I otherwise. 

then the global flux thereby defined 

= / dr n^r) (28) 

coincides with the traditional flux, through the relation II^ = n(27r/£). 

For a "good" filter obeying the condition Eq.(p5|) it follows from the (stronger) local estimate 
of the previous section that II^ = O if the minimum Holder exponent of the velocity is 



h. However, it was found in |21| that the global flux for the sharp Fourier filter may be much 



larger, and, in general, only the estimate 

n(A;) = O (k^-^^) , (29) 



may be true. In Example 1 of [^] we constructed in 3D a velocity field v G C^{Q), < h < 1, 
for which n(A;) = (const. )/c"^~^'^ at a set of /c's of arbitrarily large magnitude (namely, k = 2^ko 
with integer A > 0.) It is important here just to recall the salient features of that example. 
We have already seen that the estimate 0(2^^-^^^^^) is true if the local interactions dominate 
in the energy flux. Therefore, the example constructed was of the opposite sort, for which 
highly nonlocal triads are dominant. The idea of the construction was very simple: we used a 
Fourier series in which every octave band of wavevectors contained one conjugate pair of modes 
at the very bottom and one pair just below the top, in just such a way that the top pair of 
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each band and the bottom pair of the next highest band differ by a fixed pair of wavevectors 
in the lowest band. This pair of low-wavenumber modes induce an instantaneous transfer of 
energy between the high-wavenumber bands from the pairs at the top of one band into those 
at the bottom of the next higher band, without losing any energy themselves. This type of 
"catalytic" effect of the low-wavenumber modes for transfer between high-wavenumbers was, 
to our knowledge, first discussed in detail by Brasseur and Corssin in [Q. The same effect was 
observed in the numerical study of Domaradzki and Rogallo who described it as "local 
transfer by nonlocal triads" (although their simulation was at a rather low Reynolds number 
and did not make clear whether such interactions might persist as an inertial-range effect). The 



triads involved correspond to an "R- interaction" in the classification of Waleffe [14]. 

A remarkable feature of the example is that the fiux is proportional to the amplitude U 
of the low-wavenumber modes rather than to their gradient. It therefore corresponds to an 
interaction of strictly convective type. This interpretation of the "R-interactions" was also 



argued by Waleffe in |14], whereas Domaradzki [^] has put forward a very different picture 
"...where the role of the large scales is to produce intermittent regions of relatively strong, 
internal shears characterized by smaller length scales, which serve as regions of efficient small- 
scale transfer." The explicit verification of the convective nature of these processes renders the 
latter explanation untenable for the inertial-range at high Reynolds number. In fact, the basic 
physics of these processes was already discussed by Kraichnan in 1966 in the context of the 



energy transfer equation of ALHDIA |35]. He noted that: "Convection of high-wavenumber 
structures by strongly excited low-wavenumber velocity components implies a rapid change of 
phase of the high-wavenumber Fourier amplitudes; that is to say, a rapid exchange of energy 
between sine and cosine components of the high wavenumbers. This exchange is represented in 
(2.6) [his transfer equation] by the large, cancelling input and output contributions. The net 
contribution T<g(A;, t) represents the effect of straining alone." In other words, the convection of 
structures at a high wavenumber ~ A; by the largest, most energetic eddies will result in a rapid 
transfer of energy by a small distance in wavenumber back and forth across the sharp spectral 
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boundary at k. In a time-average of the flux (n(A;)) such effects may be expected to cancel, not 
only within the DIA closure, but in actual fact. However, for an individual ensemble realization 
at a given instant in time such effects will be present in n(A;), and, indeed, will entirely dominate 
in n(A;) at large k, swamping the comparatively small contribution of local triads. 

Although the nonlocal, convective effects may be expected to cancel in a time- or ensemble- 
average of n(A;), this would hardly suffice to justify Kolmogorov's strong statement of statistical 
independence of the small-scale modes. In fact, the previous observations have raised some 
doubts, voiced in |1C, 11, whether the small scales will really achieve such independence, 
e.g. of the external forcing or of anisotropy of the large scales. However, it seems clear from our 
discussion above that the nonlocal contributions to n(A;) are an artefact of the sharp spectral 
boundary, and do not represent an actual physical mechanism for transfer of excitation from 
large scales to small ones. We have already seen that the LES flux, with a "good" filter choice, 
gets no net contribution from such terms even at a single instant in a fixed flow realization (when 
the Holder regularity is satisfied). Instead, this fiux properly represents energy transfer as due 
to the large-scale strain. It is worthwhile to study in more detail the wavevector contributions 
to the LES fiux in order to understand the origin of cancellations in the nonlocal triads. 

The relation of the LES flux to the conventional spectral flux Il{k) is obtained very easily 
when the filter function G is spherically symmetric, as we hereafter assume. If the total resolved 
kinetic energy is represented as 

Ee{t) ^ Jdr ^v|(r,t) 
dk 1 



/ (^ ^|v(k,i)nG(fc£)p, (30) 



then using H^ = — dE^/dt , it follows that 

Euler 

= / dk U{k)De{k), (31) 
Jo 

where 

De{k) = -^\Giki)\'. (32) 
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Observe that 

/ dkD,{k) = |G(0)|2 - |G(oo)|2 
Jo 

= 1-0 = 1, (33) 

and from its definition Eq.(|32[) that 

D^{k) > 0, (34) 

if |G'(A;)p decays monotonically in the spectral radius k. Under these modest assumptions, 
therefore, Di{k) is a normahzed density in k and the LES flux II^ is a "scale-average" of the 
spectral flux Il{k). The nature of the averaging can be better seen by considering a typical 
example \G{k)\'^: 



A nice example, as pictured, will fall smoothly and monotonically from its value 1 near A; = 
to its value at k = +oo in an interval of width ^ 1/i centered at 27r/i. From this it follows 
that Di will appear as: 
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Therefore, for a reasonable choice of filter, the function will give an average over the interval 
of width ~ l/£ centered at 271 /L 

To study the nature of the cancellations involved with such an average, it is convenient to 
adopt a simple model filter G representative of this class. We take 



\G{kf 



1 



< A: < 27r 



2- 27r < < 47r , 



(35) 



47r < /c < oo 



which is pictured as: 



In this instance the averaging function is just the uniform distribution over the interval 
[27r/^,47r/£], so that 

Uf = — . dk U{k). (36) 

2vr J2TT/e 



It was this particular example which we studied in detail in |21] and it is essential to the 
arguments of this paper that the reader be familiar with that work. 



We showed in [21| that the conventional measure of spectral energy flux, n(A;), is inade- 
quate to detect the distance energy moves in wavenumber space past the cutoff wavenumber 
k. By wavenumber conservation of the triadic interactions, energy may move to at most a 
wavenumber of magnitude 2k. Efficient energy transfer of this type is achieved by the local 
triadic interactions, in which two wavenumbers of magnitude near k transfer energy to a mode 
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with wavenumber near 2k. By contrast, the Example 1 of achieved a transfer to wavenum- 
bers > A; by simply moving the energy a fixed, small distance k^ outside the sphere of radius 
k, for arbitrarily large k. The defect of the traditional flux n(A;) is that it does not distinguish 
between these two transfers, whereas a proper measure of flux should take into account the 
relative distance the energy is moved outside the sphere. Independent of any connection with 



LES, it was these ideas that led us in |21| to introduce the "band- averaged flux" 



nA = ^^^ dkU{k), (37) 



as a proper measure of energy transport. This is equivalent to the flux in Eq.(|36D if the filter 
length is chosen there as i = 2~^(27r), corresponding to a wavenumber 2^ at the bottom of the 
Ath octave band. 

By a direct argument, it was shown for this quantity in [^T[ that the "local estimate" holds 
TTa = O (2^(1-3/1)^ 

, corresponding to transfer by local triads. The main property of the Navier- 
Stokes dynamics which was used is so-called detailed energy conservation. "Partial fluxes" may 
be defined, measuring transfer from the A^th band to the Lth, induced by the Mth, as 

n^^^L = [ dr v^(r)[(vM(r)-V)v^(r)], (38) 
Jn 

in which V7v(r) is the component of the full velocity field obtained by summing only over 
wavenumber modes in the A^th octave band. Simple integration by parts shows that the partial 
fluxes satisfy 

nri. + nK^ = o, (39) 

which is the "detailed conservation" law. Because of this property, one may anticipate that 
many cancellations will occur in the expression Eq.(|3^) between the small transfers backward 
and forward between neighboring small subbands of the Ath octave band. The heuristic pic- 
ture may be explained by the following figure showing the detailed transfer of energy between 
subbands of size 2*^ for the contribution n(^^^^^(A;) to the flux induced by the Mth octave 
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band: 



This picture suggests that the "local transfers by nonlocal triads" noted by Domaradzki and 
Rogallo in |1C] will cancel instantaneously in the quantity 11^''^ ^\ The analytical verification 



of these cancellations was given in detail in Section 3 of |21|. The idea that such contribu- 
tions should cancel with an appropriate combination of triadic contributions was also pointed 



out independently by Waleffe |]1J] and Zhou |19]. What is important here is that the LES 
method leads automatically to an appropriate weighted average of the flux Il{k) for which such 
cancellations occur. 

Our previous analysis for the model filter Eq.(|35|) actually gives somewhat more than just 
the "local estimate" for = IIa. In fact, our estimates allow us to put precise bounds on the 
fraction of the flux due to distant triads. The crucial equations are Eqs.(34) & (44) of [^] which 
represent the two contributions to the flux as sums over the "deviation parameters" A, A', A" 
which measure the distance in shellnumber of the contributing shells N, M, L from A. In fact, 
it is easy to see from those equations that the fractional contribution to IIa from triads with 
smallest wavenumber in a shell lower than the (A — 5')th is 

H^glll = O (2-(i-'^)^) , (40) 

while the contribution from triads with largest wavenumber in higher than the (A + S')th shell 
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IS 



nA(5,>) 2.5V (41) 



Ha 

For the latter estimate, note that the contributions arise from terms with both L > A and 
M > A, for which A' and A" can differ at most by one because of the wavenumber selection 
rule. This gives the factor of 2 in the exponent since A' ~ A" ~ S. The parameter S measures 
the degree of nonlocalness of the triad (in IR or UV direction), i.e. the minimum number of 
"cascade steps" S between the smallest and largest wavenumbers. 

Hence, we see that distant triads in fact make an exponentially small contribution in shell- 
number to instantaneous flux when the velocity field has Holder type regularity with < /i < 1. 
This is a fact well-known for mean-flux as calculated within various analytical closures when 



the energy spectral exponent n lies in the range 1 < n < 3: see [^] and Appendix 1 |37|. 
(Note the Holder and spectral exponents are related heuristically by n = 1 + 2h.) However, 
our results here were derived for instantaneous flux in an individual flow realization without 
any closure approximation or statistical assumption. The exponent (1 — h) gives the "infrared 
locality" whereas the exponent 2h gives the "ultraviolet locality." The bound h < \ seems 
rather secure but /i > might be violated due to the appearance of "singular structures." It 
should be appreciated, however, that while the contribution of distant triads is exponentially 
small in shellnumher this corresponds only to a slow, algebraic decay in wavenumber. 

Our accounting of the fractional contribution of distant triads is very close to that made 



by Zhou in [12, 19 1. Our parameter S is essentially the same as his "scale disparity" s (more 
properly, it corresponds to log2 s) and nA(5', <), nA(5', >) are the same as his n*(A;, s), n'^(A;, s), 
respectively. In his work, however, he studied mean flux using direct numerical simulation 
at modest Reynolds number or LES to generate the raw statistics. If we take in our bounds 
Eqs.(^0|),(41) above h = ^ (the K41 mean value), then we obtain in both cases 



nA(5)/nA = o (2-2^/3j ^ (42) 

which Zhou referred to as "Obukhov scaling" in the disparity. On the other hand, he referred 
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to an estimate O ^2"'*'^/^^ as "Heisenberg scaling." In his study, both scaUngs were observed 

in certain cases, although, in our opinion, the inertial-ranges were too short to convincingly 

demonstrate either exponent. In any case, it might not be appropriate to directly compare 

our estimates with his since he dealt with an average transfer, while our flux estimates are for 

instantaneous transfer in an individual realization. Additional cancellations of distant triads 

could occur in time- or ensemble-averaging, so that better bounds on the nonlocal contributions 

could exist. What is important in our result is that, given the required Holder regularity, the 

"locality" of the energy transfer holds at each instant for a fixed turbulent flow field. 

3.2 Comparison With the ED QNM Study of Leslie- Quarini 

Superficially similar distinctions as ours above between transfer for sharp Fourier vs. graded 

filters were observed already for mean statistics within analytic closure by Leslie and Quarini 

|3§. In fact, it was observed in their work that without a low-wavenumber cutoff divergences 

appear in the sharp Fourier case for both forward transfer — "drain" — and backward transfer — 

"backscatter" — as the grid-wavenumber Ki is approached, verifying an earlier observation 

of Kraichnan based upon the TFM closure [^ ]. These divergences are due to convection 

by vanishingly low wavenumbers, and indicate a breakdown of locality for each of the two 

separate terms. However, the divergent parts cancel identically in the net transfer, leaving 

just a finite "cusp." ^ It is this cancellation of the divergences due to distant triads which 

justifies our remark above that transfer across Ki in the closures is dominated in magnitude 

by the local triads, and not by the distant triads which dominate in individual realizations for 

the sharp Fourier filter. Furthermore, for the Gaussian filter Leslie and Quarini found that 

neither forward nor backward scatter has separately the divergence (and even the finite "cusp" 

essentially disappears.) This last result is apparently in agreement with our conclusions, but a 

closer examination shows that there is no real correspondence. 

To make this point we must repeat a little of our previous analysis in the context of the 
■^In fact, they are the same divergences as those observed in Kraichnan's earher study of energy transfer in 
the ALHDIA closure cited before S. 
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closure models. We shall follow the notations of Leslie and Quarini [|3^. Within the analytical 
closures, the mean energy transfer into wavenumber k due to interaction of wavenumbers p, r 
is expressed as 

S{k\p,r) = sltk + Slt, (43) 

with 

Sl'ik = levr' • V • kVr{xy + z^)q{rMp) - q{k)]. (44) 
Here, q{k) = E{k)/A7rk'^ is the energy spectrum per Fourier mode, x,y,z are geometric factors 
associated to each triad (direction cosines) of order one, and Okpr = ['?(^) + ??(p) + ^(^)]~^ is 
an interaction time-scale associated to the triad. Different closure models are distinguished by 
their choice of r]{k), but a typical one is the EDQNM choice of r/(/c) = uk"^ + A [k^E{k))^^'^ . 
The average flux is expressed in the closures by 

roc /■ r 

U{k)= / dk' dpdr S{k'\p,r). (45) 

Jk J J Ay 

For the locality analysis, the dangerous contribution to consider is 

n«(A;)= Hdk' [ dpS^;i^, (46) 

for r <^ k, which arises from convection by the low wavenumber r. It is very easy to make a 
simple estimate of this term as 

n(^) (fc) ~ r dk' [ dp Ok'prkVrqir) [g(p) - qik')] , 
Jk JAy^ 

~ ■ k^rOrkkqir) ■ (^^) Q{k)- (47) 

The factor of arose from the fact that k',p can only take values in a range r around k. The 
crucial factor is the (r/k) which arises from expanding q{kzizr) with a power-law spectrum and 
exploiting the cancellation between input and output terms in Eg. (|4^) . 

The energy scaling corresponding to our Holder exponent h is E{k) ~ A;"^"'^^^''). This leads 
also to Orkk ~ r^~^ since the time scale for /i < 1 is clearly dominated by the lowest wavenumber. 
Hence, 
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Since the integral J^^ drr ^ is dominated by the UV region when h < 1 and is then ~ /c^ ^, it 
fohows finahy that 

n«)(A:) = drU^''\k) ~ yfci-^h^ (49) 

Jko 

exactly as expected. However, if the cancelling in-out contributions were not observed and the 
factor of (r/k) missed, then the scaling in Eq.(|48|) would be replaced by Il^'^\k) ~ /ji-2h.^-(i+h)_ 
Since the integral of r^^^^^^ is dominated by the IR region for h > 0, the false result II^^) ~ 
k^~'^^ would then be obtained. Notice here that, as is well-known, the proper locality estimates 
are obtained without any such scale-averaging as we performed previously. It would be possible 
to make such an additional scale average in the closure context, and the same cancellations will 
occur since the only requirement for our argument was detailed conservation, which is satisfied 
by the "partial transfers" : 

Sl^ik + = 0. (50) 

However, it is clear that the origin of these cancellations will be exactly the same as the "in- 
out cancellations" without averaging, since in the closures the detailed conservation condition 
Eq.(|50|) is just another consequence of the in-out cancelling terms. The scale-average is super- 
fluous for the closure flux, which is already local. Intuitively, this is clear, because the closure 
flux is supposed to correspond to an ensemble-averaged value, which will show appropriate 
cancellations between individual realizations with positive and negative values of flux. 

We can now comment on the Leslie-Quarini results on graded filters [^]. In their case. 



generalizing upon Kraichnan in |38|, they defined a wavenumber-dependent eddy viscosity 



t^nik) = VniklKi) (corresponding to Kraichnan 's i^{k\km)) as 

Un{k) = -S{k)/2k^E{k), (51) 
where they took 'E{k) = G^{k)E{k) to be the energy spectrum of the filtered velocity and 

S{k) = ^JJ^ dpdr G'^{k){l - G{p)G{r)}S{k\p,r) (52) 



21 



was what they defined as the "true subgrid transfer." Again following Kraichnan, they decom- 
posed the transfer into two terms the "drain" part, 

A 

2k''ud{k)E{k) = ^ekprbkprE{r)E{k) ■ G\k)[l - G(p)G(r)] (53) 

and the "backscatter" part 

^(^) = E ^OkprbkprE{p)E{r) ■ G^{k)[l - G{p)G{r)]. (54) 

We emphasize that these two terms arise exactly from the output and input terms of S{k\p,r), 
respectively. In particular, they do not satisfy separately a detailed conservation property such 
as Eq.(pO|), which arises precisely from a cancellation between those terms. Leslie and Quarini 
then made a corresponding division of the eddy viscosity Vnik) into Vd^k), the "drain" part, and 
Vb{k), the "backscatter" part. They numerically reproduced the result of Kraichnan that these 
two separate terms for the Fourier filter each diverges as k approaches the grid wavenumber Ki, 
and that the divergence cancels in the net viscosity Vn{k)- However, they furthermore showed 
from their numerical work that with the Gaussian filter even the separate terms, Vd{k) and 
Vbik) are finite for k approaching Ki. 

It is clear that this is what one would expect from our analysis. In fact, what we have 
established is much stronger than this, because (assuming Holder regularity) we showed that the 
local estimate holds for every individual realization with a "good" filter, whether instantaneous 
transfer is forward or backward. Since the divergences observed by Kraichnan in mean transfer 
for the sharp Fourier case are a consequence of locality breakdown, it ought to be anticipated 
that they will not occur for a "good" choice of filter. 

However, the kind of cancellations which we established for H^ in the solutions of Navier- 
Stokes cannot occur separately in the quantities Vd fmd vi, defined by Leslie-Quarini. In fact, 
there is no possibility to exploit detailed conservation of energy, as in our microscopic argument, 
because in the closure that depends precisely upon the cancellation between those terms! There- 
fore, there is no equivalence of our results and Leslie- Quarini' s. Formally, the disappearance of 
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divergences in Ud, Vb observed in |Q appear to arise in a quite different way, due to their being 
simply "smoothed out" by the averaging functions G in Eq.(^). Note that Leshe-Quarini did 
not give a theoretical explanation of their effect — as we have for the corresponding phenomenon 
in Navier-Stokes — but only observed it numerically. It seems to be, in fact, a different effect. 
The fact that our argument based upon detailed conservation does not work in the closures just 
illustrates the point that the closures do not necessarily reflect reality. In fact, it is well known 
that energy conservation does not hold realization by realization in the Langevin models of the 
closures, so that the cancellations we established due to detailed conservation in Navier-Stokes 
solutions will not occur there. 

In our opinion part of the difficulty of making a comparison of our result and that of Leslie- 
Quarini lies in the different definitions of "backscatter" in the two works. In ours here, we 
define a "backscatter" event — as in current LES practice — as a spacetime point (r, t) where 
n£(r, t) < 0. However, the Leslie-Quarini definition of "backscatter," while employing the filter 
function G, is still essentially Fourier-based. In other words, they define a "backscatter event" 
as one where, roughly speaking, the mean Fourier fiux n(A;) < and then take a suitable 
weighting with the filter function. (Even this involves an assumption that the "input" term to 
transfer arises from averaging conditioned on n(A;) < while the "output" term arises from 
averaging conditioned on n(/c) > 0: this is not obviously true.) One may expect that these 
two definitions are globally equivalent, for space-averages, but even this does not appear easy 
to demonstrate. 

4 Conclusions 

4.1 On the "Multifractal Model" 

We believe that the present work has demonstrated that the Parisi-Frisch "multifractal 
model" is not just an esoteric theoretical speculation relevant only to exotic "intermittency" 
phenomena in turbulence. In fact, it gets considerable support from the fluid equations them- 
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selves when coupled with some basic facts of turbulence phenomenology. Furthermore, the 
"multifractal model" provides a coherent framework which rationalizes and explains many ob- 
servations of experiment and simulation. Local Holder continuity is the maximal regularity 
consistent with the phenomenon of constant mean energy flux through the inertial interval. 
It is also a sufficient condition for local energy transfer, necessary for the validity of the Kol- 
mogorov universality hypothesis. While it must be regarded as tentative and so far unproved, 
we believe the "multifractal model" is valuable at least as a working hypothesis. Further inves- 
tigation, especially mathematical and theoretical, should clarify its status. 
4.2 On Locality of Transfer 

The condition that the Holder exponent lie in the range < /i < 1 is sufficient to guarantee 
locality of the instantaneous energy transfer. Nonlocal triads with wavevectors in distant bands 
from each other contribute an exponentially small fraction of the total flux. However, this 
implies only a slow, algebraic decay in Fourier space, so that very high Reynolds numbers might 
be necessary for a locally-dominated inertial-range to occur (see also 5.4 below). Furthermore, 
while local transfer is certainly necessary for statistical independence of the small scales, it 
is not sufficient. A good counterexample is Burgers equation, for which all of our analysis 
applies and in which instantaneous energy transfer is likewise local. Nevertheless, the small 
scales of Burgers turbulence do not have statistics independent of the large scales but, rather, 
largely determined by them. The key difference between Burgers and Navier-Stokes equation 
is that the former lacks the pressure force of the latter which tends to stochastically accelerate 
Lagrangian fluid particles. On the contrary, for Burgers equation the velocity of a fluid particle 
is unchanged until it is captured by a shock. The difference may be characterized as a lack of 
dynamic locality for Burgers equation, since Lagrangian time correlations will not be intrinsic 
and locally-determined in scale as they are for Navier-Stokes, but will rather be associated to 
the overall time of evolution [^]. These considerations show that some chaotic features of 
the Navier-Stokes dynamics — still insufficiently well-understood — will be necessary to produce 
independence of successive "cascade steps" and universality of the small scales. 
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4-3 On the Filtering Approach 

The filtering scheme used in the LES modelhng technique turns out to be a uniquely satis- 
factory method of scale decomposition for turbulent flows, at least in the case where a "good" 
filter is employed. (For an overview of the "filtering approach" to turbulence, see the interesting 
paper of Germano ||4^.) The requirements for a "good" filter are concrete and objective, as well 
as very easily satisfied. With such a filtering technique the basic physics of turbulent energy 
transport is made transparent. In contrast, the method of Fourier series — despite its traditional 
status in the subject — introduces many subtle and opaque effects which are just artefacts of 
using a basis of plane-waves ill-suited to describe the physical processes involved. A set of inter- 
actions like vortex-stretching which are naturally understood in physical space become difficult 
to comprehend in a wavenumber representation and are masked by other effects. 

These same criticisms apply to some apparently more sophisticated representations, such 
as wavelet bases (e.g. see [0.) Since such bases were designed to provide a simultaneous 
space-scale resolution, they might appear ideal for turbulence applications. However, we have 
found that the wavelet methods that have so far been proposed have the same problems as the 
Fourier basis and in an even more severe form! We hope to discuss the wavelet method in detail 



in a future work |42|. 
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Appendix: Multifractal Model for Navier-Stokes 

In Section 2 we reached the conclusion that Holder exponents /i < 1/3 are required to allow 
for constant mean energy flux in the inertial interval. However, it may seem to some that we 
are engaged in a mathematical "hair-splitting." After all — it might be argued — in any real 
flow there is a finite viscosity, no matter how small, and the actual velocity field, the solution 
of viscous Navier-Stokes equations, will (presumably) be smooth. However, we would like to 
emphasize that our considerations apply even for finite vq > over the long inertial-range of 
length- scales, when uq is small. It should be appreciated that the estimates we have made refer 
only to the range of scales ~ £ and they are unchanged if the velocity-field is smooth over much 
smaller length-scales. To see this clearly, consider the difference over length-scales ~ £ of a 
velocity-field with small-scale components < a removed by filtering, when a ^ ^. If v G C^, 
rj > 0, then the filtered field is C°°, or even analytic, for a nice choice of the filter function, 
and AiVa will decay rapidly in ^ for £ <C a. However, 

Aiv„(r) ^ Aiv(r), (55) 

for a. This may be demonstrated analytically using 

Aiva(r) = y ds Ga{s) [Ai+sv(r) - A3v(r)] . (56) 

Since only s < a contribute to the integral in a good approximation (exactly so, if Ga has 
compact support in the ball Ba{r) of radius a centered at r) and a <^ £, Ai_|_sv(r) Aiv(r) 

and comes out of the integral, whereas the second term is small. A quantitative estimate shows 
that 

AiVa(r) = Aiv(r) + 0(a''), (57) 

where rj may be taken to be the minimum Holder exponent of v over the ball -Ba(r). Therefore, 
the remainder is negligible if Aiv(r) ~ a ^ £, and h and rj are comparable. These arguments 
just support the intuitive idea that Aiv will not be much changed if the velocity field is smoothed 
over much smaller scales than £, whether by a mathematical operation of filtering or by physical 
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regularization due to viscosity. In a similar situation, physical Brownian paths are smooth, 
because the size of atoms a is finite, but over length-scales I ^ a they "look" like nowhere- 
difi'erentiable curves of Holder index i, which only becomes exactly true in the idealized limit 
a ^ 0. 

Our previous arguments on energy flux may be given in a form which applies to the real 
situation where Reynolds number is finite, but large. Consider a case where turbulence is 
produced in a Navier-Stokes fiuid with viscosity i^q in a box of size L by random forcing at the 
length-scales of the box-size. Suppose, in fact, that in every flow realization of the turbulent 
ensemble 



\AMr)\<voi^-j , (58) 

with some h > ^ + (5 at each point r. The constant vq has an interpretation as the typical 
variation of the velocity over the largest eddies of size ~ L. Our previous estimates then imply 
that 

„,3 / 3/1-1 

mr)\ < X (z) ' 

at each point r with probability one. However, in that case, 

Kn£)l<^flV'. (60) 



L \L, 

Let us take £u to be a scale intermediate between L and A, e.g. iu = \f\L. Then l^jL ~ (Re)""*^/^ 
in the limit vq 0. However, according to our previous discussion, (H^^,) tends to a constant 
e in the limit 0. This is not consistent with the previous estimate unless (fg) ~ (Re)^"^/^. 
We emphasize that nothing like this is observed. (At least not in three dimensions; in 2D 
simulations the level of energy fluctuations is indeed observed to rise as vq 0.) Therefore, the 
initial assumption cannot hold with 6 > for the range of scales X £ <^ L. 

The "multifractal model" accounts for the constant mean flux over this range of scales, with 
(vq) remaining finite, by postulating a suitable distribution of Holder exponents h in an interval 

[^minj ^max]) SO that 

inf [{3h-l) + i3-Dih))]=0. (61) 
27 



This is consistent with the requirement from Kolmogorov's "|-law" [^], which imphes ^3 = 1. 
The exponent for which the infimum is achieved is beheved to be ~ ^. As stated before, the 
velocity field could possibly be less regular than this, but not more so. 

These arguments imply that actual singularities will occur with some finite probability in 
the turbulent ensemble in the limit uq 0, if the currently observed finiteness of velocity fluc- 
tuations persists asymptotically to infinite Reynolds number. However, the present arguments 
do not require that singularities should apppear in finite time for Euler dynamics starting from 
smooth initial data. In fact, we consider a stationary turbulent ensemble in which the singular- 
ities may have been developed possibly over an infinite period of time. Even the consideration 
of decaying turbulence in 3D does not require true singularities in finite time. Suppose in that 
case a quasi-equilibrium state is achieved with a long inertial range of constant mean flux over 
the range of length-scales L{t) » ^ S> a(t), where L{t) is the (time-dependent) integral scale 
and a{t) is a time-dependent lower cutoff to this range. If a{t) were to decrease at a faster 
rate than the rate of decay of energy in the large length-scales ~ L{t), say, a{t) ~ e"'''*, then 
there would be little observable difference from the scenario where a{t) ^ at a finite time 
t = T^. Our previous considerations would still imply for this case that "quasi-singularities" 
must occur in the velocity field over the range of lengths L{t) ^ i ^ CL{t), to be consistent with 
the constant mean flux in that interval. 

A good example of this is the enstrophy cascade range for Navier-Stokes turbulence in 2D. 
Considerations of the same type as those above show that the enstrophy flux Z£{r) in 2D obeys 
an estimate 

Ze{r) = O [f") (62) 

in terms of the local Holder index of the vorticity field a; = V x v at point r. In the 2D case the 
natural "multifractal" picture is in fact in terms of the Holder spectrum of the vorticity field. 
Already in 1967 Kraichnan had proposed that the ultraviolet range of 2D turbulence will be an 
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"enstrophy cascade" with a constant mean flux, 

{Zi) = rj, (63) 

for r^Kr <^ £ <^ L (here L = the length-scale at which enstrophy is injected, r^Kr = the Kraichnan 
dissipation length) |^3[. From the previous estimate we see that h < will therefore be 
required. This argument is carried out in a mathematically rigorous way in a separate paper 
IQ, especially the last two parts of Section 4. The corresponding argument for the velocity 
regularity h < 1/3 of Navier-Stokes solutions in the 3D energy inertial-range is entirely the 
same, so that we may refer the reader to [44| for details. 
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